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We show a method to accelerate quantum adiabatic dynamics of wavefunctions un- 
der electro-magnetic field by developing the previous theory (Masuda & Nakamura 
2008 and 20f0). Firstly we investigate the orbital dynamics of a charged particle. 
We derive the driving field which accelerates quantum adiabatic dynamics in order 
to obtain the final adiabatic states except for the spatially uniform phase such as 
the adiabatic phase in any desired short time. Fast-forward of adiabatic squeezing 
and transport in the electro-magnetic field is exhibited. Secondly we investigate 
spin dynamics under the magnetic field, showing the fast-forward of adiabatic spin 
inversion and of adiabatic dynamics in Landau-Zener model. The connection of the 
present framework with Kato-Berry's transitionless quantum driving is elucidated 
in Appendix. 
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1. Introduction 

The technology to manipulate tiny objects is rapidly evolving, and nowadays we 
can control even individual atoms (Eigler & Schweizer 1990). And various meth- 
ods to control quantum states have been reported in Bose Einstein condensates 
(BEC) (Leggett 2001; Ketterle 2002; Leanhardt et al. 2002; Gustavson et.al. 2002), 
in quantum computing with the use of spin states (Nielsen & Chuang 2000) and in 
many other fields of applied physics. It would be important to consider the accel- 
eration of such manipulations of quantum states for manufacturing purposes and 
for innovation of technologies. In our previous paper (Masuda & Nakamura 2008), 
we proposed the acceleration of quantum dynamics with the use of the additional 
phase of wavefunctions (WFs). We can accelerate a given quantum dynamics and 
exactly obtain a target state in any desired short time, where the target state is 
defined as the final state in a given standard dynamics. This kind of acceleration is 
called fast-forward of quantum dynamics. 

One of the most important application of the theory of fast-forward is the ac- 
celeration of quantum adiabatic dynamics (Masuda & Nakamura 2010). The adi- 
abatic process occurs when the external parameter of Hamiltonian of the system 
is adiabatically changed. Quantum adiabatic theorem (Born & Fock 1928; Kato 
1950; Messiah 1962) states that, if the system is initially in an eigenstate of the 
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instantaneous Hamiltonian, it remains so during the adiabatic process (Berry 1984; 
Aharonov & Anandan 1987; Samuel & Bhandari 1988; Shapere & Wilczek 1989; 
Nakamura & Rice 1994; Bouwmeester et.al. 1996; Farhi et.al. 2001; Roland & Cerf 
2002; Sarandy & Lidar 2005). The rate of change in the parameter of Hamiltonian 
with respect to time is infinitesimal, so that it takes infinite time to reach the final 
state in the adiabatic process. However, by using our theory (Masuda & Nakamura 
2010), the target states (final adiabatic states) are available in any desired short 
time. There, infinitesimally-slow change in the adiabatic dynamics is compensated 
by the infinitely-rapid fast-forward. 

On the other hand, electro-magnetic field (EMF) is often used to control quan- 
tum states e.g., in manifestation of the quantum Hall effect (Thouless et.al. 1982) 
and manipulation of Bose-Einstein Condensates (Ketterle 2002). The acceleration 
of the adiabatic dynamics in EMF is far from being trivial, and therefore it is highly 
desirable to extend our theory of the fast- forward to systems under EMF. In this 
paper, we extend our previous theory of the fast-forward to the system in EMF, 
and derive a driving field to generate the target adiabatic state exactly besides from 
a spatially uniform phase like dynamical and adiabatic phases (Berry 1984). 

In section[5J we explain the method of the standard fast-forwarding under EMF. 
Section [3] is devoted to the fast-forward of the adiabatic dynamics with EMF. 
Examples of orbital dynamics of a charged particle are given in Section [H Section [5] 
is concerned with spin dynamics in magnetic field (MF). Summary and discussion 
are given in Section[BJ In Appendix, we shall elucidate the connection of the present 
framework with Kato-Berry's transitionless quantum driving. 



2. Standard fast-forward 

Before embarking on the fast-forward of adiabatic dynamics, we derive the driving 
electro-magnetic field (EMF) which accelerates the (non-adiabatic) standard dy- 
namics of wavefunction (WF) and drive it from an initial state to the target state 
which is defined as the final state of the standard dynamics. Hamiltonian for the 
system with electric field Eo and MF Bo corresponding to vector potential Ao is 
written as Hq = 2^(P + f-A-o) 2 - The electric and magnetic field are related to 
Ao as Eo = — and Bo = V x Ao, respectively. For simplicity of notation, 
we put - to be 1 hereafter, and we consider in the main text the case where in 
a priori there is no scalar potential V®. The results for systems with Vq is shown in 
|Appcndix A The Schrodinger equation with a nonlinearity constant cq (appearing 



in macroscopic quantum dynamics) is represented as 

.J$n ft 2 _ 2t ih ._ A . _ i hi ^ T An 

ih—^ = -- — V 2 * -^ — (V-Ao)^o Ao-VI-o + tt^o 

at Zttiq 2mo mo zmo 

-c |*o| 2 *o. (2.1) 

^o is a known function of space x and time t and is called a standard state. For 
any long time T called a standard final time, we choose ^o(t — T) as a target state 
that we are going to generate. Let ^ a (x, t) be a virtually fast- forwarded state of 
^>o ( x ) t) defined by 

!*«,(*) >= |*o («*) >, (2-2) 
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where a (> 0) is a time-independent magnification factor of the fast-forward. 

In general, the magnification factor can be time-dependent. Hereafter a is as- 
sumed to be time-dependent, a = a(t). In this case, the virtually fast-forwarded 
state is defined as, 

|*a(t) >= |*o(A(f)) >, (2.3) 

where 

A(t) = / a(t')dt'. (2.4) 



JO 

Since the generation of requires an anomalous mass reduction, we can not 
generate vI/q, as it stands (Masuda & Nakamura 2008). But we can obtain the target 
state by considering a fast- forwarded state *S?ff = ^ / f_f(x, t) which differs from fy a 
by an additional space-dependent phase, / = /(x,t), as 

*ff(*) = e v * a (i) = e^ (A(i)), (2.5) 

where / = /(x, t) is a real function of space x and time t and is called the additional 
phase. Schrodinger equation for fast-forwarded state ^ff is supposed to be given 
as 

d^FF h 2 9 ih , , ih AL 

ih—-— = V 2 * FF -- V-A FF * FF A FF • V* FF + — ^* FF 

at zttiq zttiq mo 2m,Q 

+V F f^ff - c |* FF | 2 * FF . (2.6) 

Vff = Vf F (x, t) and A FF = A FF (x, t) are called a driving scalar potential and a 
driving vector potential, respectively. Driving EMF is related with Vff and A FF 

as 

E FF = -^-^-W FF) (2.7) 
at 

B FF = V x A FF . (2.8) 

If we appropriately tune the time dependence of a (the detail will be shown later), 
the additional phase can vanish at the final time of the fast- forward T F , and we 
can obtain the exact target state 

* FF (T F ) = * (T), (2.9) 

where T F is the final time of the fast-forward defined by 

T= a(t)dt. (2.10) 
Jo 

Equation (|2.6p is rewritten with the use of equations (|2.ip . (I2.4[) and (I2.5[) as 

j_ d f -y , h 2 - T ih ,„ . , T ^A -W A 2 

-ft-j£*o + a-o — v *o-^ — (V-A )1<o - + 77^*0] 

at 2mo zrao mo 2mo 

[V 2 * + «(V 2 /)* + 2zV/ • W - (V/) 2 * ] 



2mo 

ih . s T ifi . r .,_„ iT „ T -, A 2 



-(V • A FF )* A FF • {i(V/)* + V^o} + tt^^q 

2mo mo 2mo 

+I/ FF *o + (a-l)c |* | 2 *o, (2.11) 
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where /(x, t), \E , (x, A(t)), a(t), A (x, A(t)), Ayjr(x, t) and Vff(x, i) are abbre- 
viated by /, c*, Ao, Afp and Vff, respectively. The same abbreviation will 
be used throughout in this section. Real and imaginary parts of equation (|2.11[) 
divided by ^/q yield a pair of equations as 



2 v. ( v/- ftA °: A " ) + 2R e [^ v^](v/- aAo : A ^ ) 

n n 

+ (a- l)Im[* V 2 *5] = 0, (2.12) 



and 



77io » 2mo 2mo at/ 

'' A " V/-(a-l)?|* | 2 . (2.13) 



777-0 ft ^ 

We can take the driving scalar potential from equation (|2.13p and the additional 
phase / which is a solution of equation (|2.12[) . 

(a) Additional phase and driving field 

In order to obtain the driving field, we should first calculate the additional phase 
in equation (|2.12p . Here we derive a general solution of equation (|2.12p from the 
continuity equation for ^o and ^ff- With the use of equation (|2.1j) . we have a 
continuity equation for ^o 

|l*ol 2 = ■ (M* v*5] X |VJ, ° |2) ' (2 - 14) 

and by using equations (|2.5p and (|2.6p . the continuity equation for ^ff is 

^|*ff| 2 = —V • (-V/|* | 2 + Im[* Wo1 - ^l^ol 2 ), (2.15) 
at mo n 

where ^ff(^, t) is abbreviated by ^>ff- From equation (I2.5P which is the definition 
of we have a relation between time derivatives of V^o and ^ff as 

||*FF| 2 =a||*o| 2 . (2.16) 

Combining equations (|2.14[) . (|2.15l) and (|2.16[) . we have the gradient of the addi- 
tional phase 

V/(x, t) = (a- l)Im[V* /*o] + (a 4r " ^F)- ( 2 ' 17 ) 

Noting the equivalence of gauges for Ao and Aff due to the initial condition 
(A (< = 0) = Aff(£ = 0)), we can take the gradient of the additional phase and 
Aff as 

V/ = (a-l)Im[Wo/*o], (2-18) 
Aff = aA . (2.19) 
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And it is easily confirmed that equations (|2.18p and (|2.19[) satisfy equation (|2.12p . 
Equation (|2.19[) implies that MF should be magnified by a times, that is, 

B FF (x, t) = a(t)B (x, A(f)), (2.20) 

where Bo = Vx Ao and Bj?f = V x A.pp. When the standard state is written 
with its phase ?7(x, t) as 

*„(x,t) = |*o(x,t)|e i, '<*'*>, (2.21) 

equation (|2. 18[) leads the expression of the additional phase and its gradient as 

V/(x,t) = (a-l)Vj7(x,A(t)), (2.22a) 
/(x,t) = (a-l)»j(x,A(t)). (2.22b) 

In equation (|2.22bj) . a space- independent constant term was neglected. 

Substitution of equation (|2. into equation (|2.13[) yields the driving scalar 
potential 

~—Vf - Im[V*„/*„] - -^(V/)> - - V/ 

mo 2mo mo n 

where / and V/ are given by equation (12.22[) . Therefore, once we have and Ao, 
the driving field can be obtained from equations (|2.7[) and f|2 . 8[) with the use of 
equation f|2 . 19[) and (|2.23[) . By applying Epp and B FF against the initial standard 
state, we can generate the target state in any short time Tp related to the standard 
final time through a by equation (|2.10[) . 



3. Fast-forward of adiabatic dynamics 

So far, we presented the fast-forward of the standard dynamics in EMF which 
enables to generate the target state in any desired short time. Here we show the 
fast-forward of adiabatic dynamics of WF under EMF in an analogous manner 
which we employed in the previous paper (Masuda & Nakamura 2010). In the 
process, the regularized adiabatic dynamics is introduced as standard states and we 
fast-forward them with infinitely large magnification factor a. In the acceleration, 
the initial state is stationary and it becomes back to stationary state at the end 
of the fast-forward. However we can not directly apply the theory in last section 
to adiabatic dynamics, because the stationary state is just an energy eigenstate of 
the instantaneous Hamiltonian and, therefore, is not suitable to be fast-forwarded. 
Thus we first need to regularize the adiabatic dynamics for the fast-forward. 

Let us consider vf^ under Eo and Bo corresponding to the vector potential 
Ao = Ao(x, R(t)) which adiabatically varies, where R = R(t) is a parameter which 
is changed from constant Ro as 

R(t) = R + et. (3.1) 
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The constant value e is the rate of adiabatic change in R(t) with respect to time 
and is infinitesimal, that is, 

dR(t) . . 

~d~T = ' (3 - 2) 

e < 1. (3.3) 

The EMF is changed through this parameter. Hamiltonian of the system is repre- 
sented as 

(p + Aofx,^))) 2 

H o = 7, 1 (3.4) 

ZrriQ 

and Schrodinger equation for is given as 

where Co is a nonlinearity constant. The results for systems with Vq are shown in 
|Appcndix A| If a system is in the n-th energy eigenstate at the initial time, the 
adiabatic theorem guarantees that, in the limit e — > 0, \&o remains in the n-th 
energy eigenstate of the instantaneous Hamiltonian. Then Vl^o is written as 

* (x,t,i2(t)) = M^R(t))^ I ° En(Rlt ' ))dt 'e im , (3-6) 

where E n = E n (R) and <f> n — n (x, R) are the n-th energy eigenvalue and eigenstate 
corresponding to the parameter R, respectively, and T = T(t) is the adiabatic phase 
given by 

r(t) = *y J dxdtfa-fa, (3.7) 

which is independent of space coordinates. <p n fulfills 

t = »• ^ 

ti 2 ih ih A 2 

-V 2 0„-- — (V-A o )0„ - — A - V0 n + —5-^) n -co|0„|Vn = ^ ^ 



2mo 2mo m.o 2mo 



(3.9) 



The second and the third factors of the right hand side of equation (|3.6|l are 
called dynamical and adiabatic phase factors, respectively, which are also space- 
independent. (We will not intend to realize these phase in the fast-forward.) The 
adiabatic dynamics in the limit e — > takes infinitely long time until we obtain an 
aimed adiabatic state (target state). 

For the fast-forward of the adiabatic dynamics, we should first choose an ap- 
propriate standard state and Hamiltonian. The original adiabatic state is not ap- 
propriate as the standard state. Quantum dynamics in equation (|3.5[) with small 
but finite e inevitably induces non-adiabatic transition, but tyg in equation (|3.6p 
ignores such transition. To overcome this difficulty, we regularize the standard state 
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and Hamiltonian corresponding to the adiabatic dynamics (Masuda & Nakamura 
2010), so that the following two conditions are satisfied. 

1. A regularized standard Hamiltonian and state of the fast-forward should agree 
with Hq and except for space- independent phase, respectively, in the limit e — > 0; 

2. The regularized standard state should satisfy the Schrddinger equation corre- 
sponding to the regularized standard Hamiltonian up to 0(e) with small but finite 
e. 

Hereafter ^f^ 9 ^ and H^ 9 ^ denote the regularized standard state and Hamilto- 
nian, respectively, which fulfill the conditions 1 and 2. In the fast-forward, we take 
the limit e — > 0, a — > oo and ae ~ 1. Applying this regularization procedure in 
advance, the adiabatic dynamics 4> n (R(0)) — > (f> n (R(T)) can be accelerated and the 
target state <p n (R(T)) is realized in any desired short time, where T is a standard 
final time which is taken to be 0(1/ e). 

(a) Regularization of standard state 

Let us regularize the standard state so that it can be fast-forwarded with in- 
finitely large magnification factor. Let us consider a regularized Hamiltonian Hq T 



H ° ~ 2^ + V ° ' (3 - 10) 



The scalar potential Vg res ^ in the regularized Hamiltonian is given as 

vt 9 \^t)=eV(^t). (3.11) 

V is a real function of x and t to be determined a posteriori, which is introduced 
to incorporate the effect of non-adiabatic transitions. On the other hand, for the 
vector potential we put A( reff )(x, t) — Aq(x, i). It is obvious that H^ 9 ^ agrees 
with Hq in the limit e — > 0, that is, 

lim# (res) (x,f) = H (x,R(t)). (3.12) 

The standard state in the adiabatic dynamics should fulfill Schrddinger equation 
up to 0(e). We suppose that a regularized standard state is given by 

^) = ^ e -i/ *B n (fl(O)*' e tefl j (3 13) 

where 9 = 0(x,t) is real, and 4> n — (f) n (x,R(t)) and E n — E n (R(t)) arc the n- 
th energy eigenstate without dynamical phase factor and eigenvalue of the original 
Hamiltonian Hq. cj> n satisfies the instantaneous eigenvalue problem in equation (|3. 91) . 
The Schrddinger equation for regularized standard system is represented as 

= _-^v 2 < e9) --^(V-A )^ re5) -— A • V< e9) + #2_< e9) 
at 2mo 2m m 2m 

+eV^ re9) - c |< e9) | 2 < e9) . (3.14) 

Substituting equation (|3.13|) into equation (|3.14l) and eliminating the equation of 
0(1) with the use of equation p.9[) . we find the equation for 0(e): 

- h ~%^ = -17— i 2 ^ 9 ■ V <A" + *(V 2 #)0„] + V4> n + — A • (V0)0„. 
OK at Zitiq Too 

(3.15) 
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Multiplying equation (|3.15[) by 10* and taking the the real and imaginary parts of 
the resultant equation, we have 

\<j> n \ 2 V 2 9 + 2Re[CV0„] • W + ^Re[c||] = 0, (3.16) 

\ = -^m'M ~ ^<t>n/<t> n ] ■ V* - ■ VO]. (3.17) 

From equation (|3.16p . 9 turns out to be dependent on t only through R(t). Therefore 
the minor term d9/dt(= e-jm) was suppressed in Ea (|3.17p . Equations (|3.16p and 
(|3.17p give 9 and V, respectively. It is worth noting that 9 is not explicitly affected 
by EMF. 



(6) Additional phase and driving field for fast-forward of adiabatic dynamics 
The regularized standard state in equation (|3.13[) is now written as 

j# e »> = \0 n \ e ^ +e9) e-* /o E ~ d \ (3.18) 

where r\ is defined as a phase of <p n by <p n — \(f) n \e tri . By using ^f^ 9 ^ in equation 
(13. 18)) instead of in equation (|2.12[) , we have 

l^| 2 V • (V/ - ^L-AlL) + 2[|0„|V|0„|](V / - ^L-AlL) 

-(a - l)[2|0 n |V(r? + e9) • V|0„| + |^| 2 V 2 (r/ + e9)] = 0, (3.19) 

where 0„(x, R(A(t))), /(x, t), A (x, R(A(t))), A FF (x, t), r?(x, i?(A(i))) and 0(x, ii(A(t))) 
are abbreviated by <p n , /, Ao, A.pp > V an d $ respectively, and the same abbrevia- 
tions will be taken hereafter in this section. 

Multiplying </>* on both sides of equation (|3.9p and taking its imaginary part 
with the use of <j> n — |</> n |e 117 , we have 

— [2|0„|Vt7 • V|^| + V 2 7# n | 2 ] + A[ V . ^|0„| 2 + 2|0„|V|0„| ~] = 0. 
mo mo n n 



(3.20) 



With the use of equation (|3.20p in equation (|3.19p , we obtain 

i ,— , A - Ajf F n/TTjr A -A FF . 



knl'V • (V/ - u fc + 2[K|V|^|](V/ - 
We can easily confirm that 



(a - l)e[2|6,|V0 • V|0„| + |0„| 2 V 2 0] = 0. (3.21) 



V/ - A ° J^ FF = (a - l)eVfl (3.22) 
h 

satisfies equation (13.211) . Noting App(t = 0) = Ao(t = 0), we have the vector 
potential App and gradient of the additional phase from equation Q3.22p as 

A FF (i)=A (A(i)) (3.23) 
V/ = (a - l)eVd, (3.24) 
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which should be compared with the result in equations (I3.8[) and (|3.9p in the case 
of the standard fast-forward. It is noteworthy that we do not have to magnify the 
MF for the fast-forward, while in the standard fast-forward we need to magnify the 
MF by a times as shown in section [2] The result in equation (|3.22[) is also obtained 
from the continuity equation. 

As mentioned in Section [3J in the standard fast-forward with a standard scalar 
potential Vq, the driving scalar potential is given by equation ([A ip . By using 
in equation (|3.13[) and V^ re9 ^ in equation p. lip instead of and Vq and noting 
equations (|3.9p . (|3.17p . (|3.23p and (|3.24p , equation ([A ip leads to the driving scalar 
potential as 

Vff , * \ En da „ , , 90 h , , ,_ „. , h An „ . 

-aelm[^/<f> n ] - as —Im[^] ■ V6, (3.25) 
on m (p n 

where we omitted a term of O(e). While the first term diverges with infinitely large 
a, it concerns only with spatially uniform phase of WF, which we do not care about 
in the fast-forward and can be omitted. Consequently, we have the driving scalar 
potential 

Vff da . , , 89 h 9 h A 

= — -e6 - qTe 2 — - - — a 2 e 2 (V6» 2 as— - ■ V6» 

n at on 2mo mo n 

-OElm^w/M - <*E — Im[^] • V0. (3.26) 

on m (p n 

The driving field can be obtained from equations (|2.8p , p.23p and p.26p . So far we 
considered the fast-forward in the systems without scalar potential Vq . The driving 
scalar potential for systems with Vq is shown in |Appendix A while the driving 



vector potential has the same form as in systems without Vq. 

The present theory of the fast-forward is different from the reverse engineering 
approach based on the inverse technique {e.g., Palao et.al. 1998) in the following 
sense: the latter approach is concerned with a direct fast-forward of the adiabatic 
state itself. By contrast, we consider the adiabatic states except for the spatially uni- 
form phase but together with a controllable additional phase, and combine an idea of 
the infinitely-fast-forward and infinitesimally-slow adiabatic dynamics. Muga et.al. 
employed the reverse engineering approach to accelerate the adiabatic squeezing 
or expanding of BEC wave packet (Muga et.al. 2009; Chen et.al. 2010) in a tun- 
able harmonic trap, finding a promising time dependence of the trapping frequency. 
While their approach is limited to WF under the harmonic trap, the present theory 
of the fast- forward enables to accelerate adiabatic dynamics of WF in any potential 
and EMF. 



4. Examples 

So far we showed theoretical framework of the fast-forward of adiabatic dynamics 
in EMF. Here we give some examples of the fast-forward of adiabatic dynamics 
without nonlinearity constant (cq = 0) for simplicity. Our purpose is the realization 
of target states defined in the adiabatic process in any desired short time, while 
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the target states are reached through infinitely long time in the original adiabatic 
dynamics. In the following examples, the magnification factor is commonly chosen 
(for < t < Tp) in the form 

a(t)e = v(l - cos(^t)), (4.1) 

where v is time average of a(t)e during the fast-forwarding, and the final time of 
the fast-forward Tp is related to the standard final time T with v as Tp — eT/v (see 
equation (|2.10p ). eT and Tp are taken as any finite value, although e is infinitesimal 
and T is infinitely-large. Namely we aim to generate the target state in finite time, 
while the state is supposed to be obtained after infinitely long time T in the original 
adiabatic dynamics, ae starts from zero and becomes back to zero at the end of 
the fast-forward. 



(a) Fast-forward of adiabatic squeezing of wave packet with electro-magnetic field 

We consider an adiabatically squeezed wave packet (WP) in two dimensions 
under the adiabatically increasing MF 

B = (0,0,R(t)), (4.2) 

where R(t) = Ro+et as given in equation (|3.ip . The vector potential corresponding 
to the MF can be taken as 

A = (-R(t)y, 0,0). (4.3) 



The lowest energy eigenstate (j) n =o of the instantaneous Hamiltonian with energy 
hr 

2mo 



jR(t) 2, 2, H(t).» 

(j) = 4/_LZ e 4R ( X +V > e l 2K . (4.4) 

V 2nn 

The corresponding regularized standard state is written with the phase of O(e) as 



This WF is squeezed when MF is increased. We accelerate the manipulation which 
control the width of WP. Equation p,16[) for is rewritten as 



_■, ; R d6 R d0 mo 2 2\ m o 



It can be easily confirmed that 



AhR 



(x 2 + y 2 ) (4.7) 



satisfies equation (14.61) . Equation (13.17[) leads to the regularized standard potential 
as 

V^ re9) =sV= -e^-. (4.8) 
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With the use of equations (|4.5I) and (14.71) in equation (I3.26|) . the driving scalar 
potential is obtained as 

Vff ,da rrin 9 9 3mn ,,99, xy , , . 

— = - a£ ^ x + y) ~ °*£> (4 - 9) 

where we omitted a spatially uniform term because it is concerns only with the 
spatially uniform phase. The driving EMF can be obtained from equations (|2.7[) 
and (|2.8p with the use of equations (|3.23l) , (I4.9[) and (|4.3[) . For numerical calculation 
the parameters are chosen as ^ = 1.0, T F = 1.0, v = 1.0, R = 1.0. The WP profile 
is shown in figure [T] at the initial (upper figure) and final (lower figure) time of the 
fast-forward. It can be seen the WP is squeezed successfully. To check the accuracy 




Figure 1. Wavefunction |$ff| profile before (upper figure) and after (lower figure) the 

squeezing. 

of the acceleration, we evaluated the fidelity which is defined by 

F = I < * FF (t)|* (A(t)) > I, (4.10) 

i.e., the overlap between the fast-forwarded state \l/ FF (i) and the corresponding 
standard one ^>o(A(t)). It is unity when ^pp(t) — ^o(A(t)). We confirmed that the 
fidelity first decreases from unity due to the additional phase / of the fast-forwarded 
state, but at the final time it becomes unity again (see figured]), which means the 
exact fast-forward of the adiabatic state besides from the spatially uniform phase 
factor. 
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Figure 2. Time dependence of fidelity. 

(6) Fast-forward of adiabatic transport in electro-magnetic field 

Here we show the fast-forward of adiabatic transport of WF in 2 dimensions 
subjected to EMF, without leaving any disturbance on the WF at the end of the 
transport. 

WF takes a form ijj(x)e~Ti E7lt , which is stationary except for the adiabatic phase, 
in the presence of vector potential A(x, y, z) at the initial time. The MF is adiabati- 
cally shifted with infinitesimal velocity e in x-direction. The shifted vector potential 
Ao (x, t) are represented with the use of A as 

A = A(x-et,y,z). (4.11) 

The corresponding regularized WF is supposed to be given in the form as 



(reg) 



(4.12) 



<f> n is the n-th energy eigenstate of the instantaneous Hamiltonian written as 

4> n = ip(x - et,y,z) = ip(x - R(t),y,z). (4-13) 

R(t) which characterizes the position of WP in x direction is adiabatically changed 
as R(t) — Rq + et with e <C 1 and Rq = 0. Equation (|4.13[) leads to the relation 



dR 



dx 



(4.14) 



In the same manner as used in the previous example, we can obtain 6 and V as 

n TO 



V = -A, 



Therefore the regularized potential V^ 9 ^ (x, y, z, t) in equation p. lip is given as 



(4.15) 



(reg) 



-eA x 



(4.16) 
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With the use of equations (I4.14j) and (14. 15)) in equation (|3.26p , we have the driving 
scalar potential 

dct 

V FF (x, y, z,t) = -m —ex - aeA x , (4-17) 

where we omitted the spatially uniform term because it concerns only on the spa- 
tially uniform phase of WF. The driving MF is shifted with time and the corre- 
sponding vector potential is given from equation (13.23[) as 

A FF (x,y,z,t) = A(x-R(A(t)),y,z). (4.18) 

The driving electric field is given by 

E FF (x,y,z,t) = ea(-^)-VV FF . (4.19) 

It is worth noting that since we have derived 9 without giving any specific profile 
on <fi n , the formulas of the driving scalar potential in equation (14. 17)) and driving 
vector potential in equation (14. 18)) are independent of the profile of the WF that 
we are going to transport. The resultant electric field due to this term and the time 
derivative of A F F can be interpreted as balancing with Lorenz force perpendicular 
to the transport. 

As a concrete example of such accelerated transport, we consider a case that a 
WP trapped in uniform MF B = (0, 0, B). We choose the vector potential as 

A = (0,B(x-et),0), (4.20) 

which leads to Bo and electric field of 0(e) in y-direction: 

E = -^ = (0, £j B,0). (4.21) 

In this case, V^ re9 ' = as seen from equation (|4.16[) . The WP is adiabatically 
moved due to 9 and Eo, while the MF does not change. From equations (|4.18[) and 
(|4.20j) . it is obvious that we do not have to change MF for the fast- forward. An 
eigenstate with energy Eq — of the instantaneous Hamiltonian with R is given 



<fi n {x,y,t)e l2m », (4.22) 



with 



^ J^e-^-mr + v\-i^^. ( 4. 2 3) 
V 2nn 

Note that <fi n in equation (|4.23p is a stationary state with an instantaneous value 
of R. We transport this state by the driving field. In this case, the driving scalar 
potential in equation (|4. 17)) is represented as 

dot 

V FF {x, y, z, t) = -m —ex, (4.24) 
at 
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which leads to the driving electric field in x-direction. The driving electric field in 
y-direction given by 

E M = _^F =a£B (4 25) 

is balancing with Lorenz force in classical picture. Without this term the path of 
WP would be bent in y-direction. 

In the numerical calculation the parameter are chosen as ^ = 1.0, Tp = 1.0, 
v = 8.0 and Ro = 0. By applying the driving potential in equation (|4.24p . we 
accelerate WP. In figure [3J the WP profile |\&ff| is shown at the initial and final 
time of the fast-forward. The WP is transported by distance 8.0 in time 1.0 and 
becomes stationary at the end. We confirmed that WP is moved without changing 
its amplitude profile during the acceleration. We evaluated the fidelity defined by 




Figure 3. Wave packet profile [\PffI at initial and final time of the fast-forward. 

equation (|4.10l) and confirmed that it becomes back to unity at the end of the fast- 
forward (see figure 2]) . Thus we have obtained the adiabatically accessible target 
state in a finite time Tp — 1.0. The inset in figure 0] shows the path which the WP 
traces under the correct driving field (thick line) and the path without electric field 
in y-direction aeB (thin line). Without this electric field, the orbit of WP is bent 
in y-direction by MF. 

Here we considered the case that WP is trapped by MF. In |Appendix A| the 
driving field is shown in the case that WP is trapped by both MF and a scalar 
potential Vq. 
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Figure 4. Time dependence of fidelity. The inset represents the path which the wave 
packet traces under the correct driving field (thick line) and one without electric field in 
y-direction (thin line). 

5. Fast-forward of spin dynamics in magnetic field 

We now proceed to the dynamics of spin with S = \, by suppressing the orbital 
degree of freedom. The fast-forward of spin dynamics in MF is as follows. Let us 
consider the Hamiltonian given by 



Ho(t) = \{Bf\t)a x + B?\t)a v + (*)*,), 



(5.1) 



where B^ , By'' and B^'' denote components of the MF and a x , a y , a z are Pauri 
matrices. In equation (|5.1[) the negative of Bohr magneton was suppressed for sim- 
plicity. The MF can change with time t. Schrodinger equation can be given in 
<r z -diagonal representation as 




„(0) 

i 

„(0) 



where 



H lx {t) 
H 12 (t) 



\Bf\t) 

\{B^\t)-iB^{t)). 



The state 



*o(t) 



.,«>) 

,(0) 



(5.2) 

(5.3a) 
(5.3b) 

(5.4) 



is defined as a standard state that we are going to fast-forward. The exactly fast- 
forwarded state ^ a (t) is defined by \&o(A(i)) with A(t) in equation (12.41) . Let the 
Hamiltonian for \I/ a to be represented as 



H a = 



H ll 



2 J 




(5.5) 
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The definition of the exactly fast-forwarded state in equation (I2.3|) leads to the 
relation 



<N> r 



dt 



at t=A(t') 



With the use of Schrodinger equation for ^> a with Hamiltonian H a together with 
equations (15.2[) and (15.6[) . we have the relation between H a and Hq as 



H a (t)=a(t)H (A(t)). (5.7) 
Therefore we can obtain the driving MF, Hff, as 

B FF (t)=[ By F (t) = «(t)i?< 0) (A(*)) • (5.8) 
V Bf F (*) / V a(t)B^(A(t)) j 

In the fast-forward of spatially distributing WF, there was a problem of the 
anomalous mass reduction. To resolve this problem we introduced an additional 
phase / on the fast-forwarded state (see equation (|2.5p ). but here in spin dynamics 
we do not have to use it. And we just require to magnify the MF to generate the 
exactly fast- forwarded state 

(a) Fast forward of adiabatic spin dynamics 

Now, we show the fast-forward of adiabatic spin dynamics. Suppose that Hamil- 
tonian which is adiabatically changed is represented as 

H{R{t)) { HUR{t)) -H n (R(t)) ) ' (5 ' 9) 

where R(t) is the parameter which is adiabatically changed as in equation (|3.1[) . 
namely, e<1. The matrix elements of the Hamiltonian are related to those of MF 
as. 

Hu(R(t)) = l -Bf\R{t)) (5.10a) 
H l2 (R{t)) = l(Bi°\R(t))-iB^(R{t))). (5.10b) 



Suppose 



* (i*(*)) = ( llf^] ) e -*JS-Wt))* 6 ir(t) i (5 . n) 



to be an adiabatically evolving state. T(t) is an adiabatic phase, which is common 
to both component. We have 

e <«» ( :S«! ) - ( ) - ^ 

for an instantaneous Hamiltonian with a parameter R. In this fast-forward, we can 
not utilize the same manner as used in the standard fast-forward of spin dynamics 
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because the MF diverges due to infinitely large a (— 0(l/e)). Thus we need to 
regularize the standard Hamiltonian. 

As we did in the fast-forward of adiabatic orbital dynamics, we regularize the 
system so that the Schrodinger equation is fulfilled up to 0(e) for regularized stan- 
dard state and Hamiltonian. In this case, however, there is no spatial distribution 
of WF, and therefore we do not have to regularize WF with any additional phase of 
0(e) corresponding to 9 in equation p,13[) . We simply put a regularized standard 
state without adiabatic phase T as 

*^ } W*))=(^jj )e-*J?*W*»* (5.13) 

for which the regularized Hamiltonian is given by 

n {reg) m(A s_ f H n (R(t)) + sh n (R(t)) H 12 (R(t)) + eh 12 (R(t)) \ 
{ [>> \ H* 12 (R(t)) + eh* l2 (R(t)) - Hll (R(t)) - eh lx (R(t)) J' { ' 

introducing additional terms h\\ and h\ 2 . Schrodinger equation is written as 



ih^^-=H$ ea) *§ e ° ) . (5.15) 



The use of equation (|5.12p in equation (|5,15[) leads 



3c\ 

ift-Qjj = huci + hi 2 c 2 , (5.16a) 
dc 2 

ih— = h* 12 ci - hnc 2 . (5.16b) 



dR 

From equations (|5 . 1 £>[) we can obtain 

hn = ih{c " l JS + C2 l^ ) - (5 - 17a) 

where h\\ in equation (|5.17a[) is purely real because ^(|ci| 2 + | c 2 1 2 ) = 0. 
The fast-forwarded state is given as 



FF 



t Cl (R(A(t))) \ ± ft E (R(A(t)))dt (5 lg N 

W { c 2 (R(A(t))) ) ■ [ ° A *> 



Since there is no additional phase on f j?f , it is the exactly fast-forwarded state 
except for the adiabatic phase common to both components. We suppose that 
Hamiltonian Hpp drives ^ ff- The time derivative of equation (|5.18[) is given by 

where *f.f(*)> a(i), E(R(A(t))), c x (i?(A(i))) and c 2 (R(A(t))) are abbreviated by 
a, E, ci and c 2 , respectively, and the same abbreviations are used hereafter 
in this section. Schrodinger equation: 

ih ^dT~ = Hff ^ ff ( 5 - 2 °) 
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and equations (|5.12l) . (I5.16[) and (|5.19p lead to the driving Hamiltonian as 



HfF (t) - r U {FF)\ , U {FF) U {FF) 




B iFF) B (FF) ^iFF) 



iBy ' -B\ 

H12 + aehi2 
— if 11 — aehu 



(5.21) 



where Hu(R(A{t))), H 12 (R(A(t))), h n (R(A(t))) and h 12 (R(A(t))) are abbreviated 
by if 11, ifi2, hn and h\2, respectively. The driving MF is obtained from equations 
(jSTTU)) . (jCTfl) and (15T2"Tj) as 

B (FF) = Bi o ) _ 2hEalm{c *^±_ Cl ^l )) (5 . 22a) 

< F) = B r- 2feaRe ^M- Cl tf)' (5 ' 22b) 
fl(") = s(°)-2fe a Im( C ^ + C2 ||) ) (5.22c) 

where B\ (t) and B\ (R(A(t))) are abbreviated by B\ and B\ , respec- 
tively, and i denotes x,y,z. Noting set — 0(1), the excess field, i.e., the differ- 
ence between Bj?f and Bo gives a nontrivial contribution. By applying Hff = 

(pp) ( F F*] ( F F*] 

(£?i , -By , -Bz ; ) given in equation (15.221) . we can accelerate adiabatic dy- 
namics and obtain, in any short time, the exact target state except for the common 
phase between both components. It is obvious that ^>ff and Hff coincide with 
ancl Bo, respectively, at the initial and final time. 

In closing this subsection, it should be emphasized: Our purpose in this paper 
lies in the fast-forward of the adiabatic dynamics except for the uniform phase. 
Therefore we suppressed the adiabatic phase in equation (|5.13l) . Equation (|5.22[) is 
a result of such a simplified procedure. On the other hand, we can also accelerate the 
adiabatic dynamics with the adiabatic phase being included in equation (j5.13[) . The 
resultant driving field is slightly different from equation (|5.22l) . which is described 
in |Appcudix B| There, the driving field proved to be equal to the one obtained 
recently by Berry (Berry 2009). 



(6) Examples 

As an example we consider the adiabatic dynamics in which MF is rotated 
adiabatically into opposite direction while its magnitude is kept constant. Let the 
MF be written as 

(sin9(t) cos ip \ 
sin0(i)sin</J , (5.23) 
cos (9(f) / 

where 

6{t) = R(t). (5.24) 
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R(t) is given by equation (|3.1I) with Rq = 0. B and (p remain constant. e r is a unit 
vector pointing the direction of MF. Hamiltonian for this MF is represented as 

H(R( ^ B ( coa(R(t)) sm(R{t))e-*r \ , . 

- 2\sm(R(t))e^ - cos(R(t)) J' ^ j 



Under this magnetic field, the regularized standard state corresponding to an adi- 

B 
2 



abatic state with eigenvalue A + = is given by 




< r£9) = ( ~". \t) (5.26) 
\ e v sin — j 2 - 

From equations (|5.17|) . (I5.24p and (|5.26[) . we have 

/in = 0, (5.27a) 
h\2 = — g ( sm V + * cos V) ■ (5.27b) 
From equations (JEHJ), (f5T22|) . and ([5~27)l . we can obtain the driving MF as 

B i FF) (t) = B sin{R(A(t))) cos (p - £a(t)H siniy9, (5.28a) 
B { v FF) {t) = B sin(R(A(t))) sin (p + ea(t)H cos (p, (5.28b) 
Bi FF \t) = Bcos(R(A{t))), (5.28c) 

where A(t) is given by equation (|2.4I) . In numerical calculation the parameters are 
chosen as v = tt, T f = 1.0 ^ = 0.0, R = 0.0, fl(T F ) = ir and B = 1.0. Therefore in 
the fast-forward, the direction of MF switched into opposite direction in time Tp. 
Initially the state is set as 

*FF(t = 0) = ( J V (5.29) 
Spin is rotated by the applied B FF and points the opposite direction (?) at the 



v 1 t 

end of the fast-forward. In figure [51 the dynamics of the spin state: |ci Fi? ' > (t)\ 2 and 

k^^WI 2 are shown. It is confirmed that spin is flipped exactly from up to down 
at the final time Tp- The spin state becomes stationary after the acceleration again. 
We also confirmed that the fidelity defined between ^pp(t) and *&o(A(t)) is unity 
during the fast-forward because there is no additional phase on the fast-forwarded 
state. We can control the direction of spin more generally, in any desired short time, 
by making tp changing through R like 9. 

As another example, we show the fast-forward of adiabatic dynamics in Landau- 
Zener (LZ) model (Landau 1932; Zcncr 1932). We consider MF: 



B(i) = , (5.30) 
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Figure 5. The time dependence of Ic^^l 2 (solid line) and Ic^^l 2 (broken line). 

where A is a constant, and R(t) is given in equation (|3.ip with large negative 
constant i?o- The Hamiltonian is given by 



2 V A -R(t) 
The adiabatic state with eigenvalue A+ = ^ R2 + A ' 2 is given by 

*t 9 \t)=( R -^ )e-i^, (5.32) 

where 



Here we have 



s = {2^/R 2 + A 2 (\/i? 2 + A 2 - R)} 1 ' 2 . (5.33) 
da 1 A 



OR 2V2Q 5 / 2 
dc 2 1 (Q- R) l ' 2 {Q + R) 



OR 2V2 Q 5/2 



(Q-R)* (5.34a) 

(5.34b) 



where 



Q = yjR 2 + A 2 . (5.35) 

By using equations (15.32)) and (|5.34l) in equation (|5.22l) . we have the driving field 

as 

B FF (t) = ( -ea{t)h R(Ht f )2+A , \ , (5.36) 

V R(Mt)) J 

The y-component of Hp F is identical to the one appearing in Landau-Majorana- 
Zener model (Landau 1932; Majorana 1932; Zener 1932; Berry 2009). The dynamics 
of the spin state driven by T5 FF is shown in figure [6] during the fast- forward. The 
parameters were chosen as v = 100.0, T F = 1.0, Rq — —50.0 and A = 1.0. We 
confirmed the fidelity is unity throughout the acceleration. Interestingly, with the 
use of the formula different from equation (|5.22p . Berry also obtained equation 
(I5.36[) . This mystery is solved in |Appcndix B| 
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Figure 6. The time dependence of |cj 



(PF)|2 



(solid line) and |c 



(broken line). 



6. Conclusion 



We have presented the theory of the fast-forward of quantum adiabatic dynamics 
in electro- magnetic field (EMF) . We derived the driving EMF which accelerates the 
adiabatic dynamics and enables to obtain the final adiabatic states besides from the 
spatially uniform phase in any desired short time, while the final state is accessible 
after infinite time in the adiabatic dynamics. In the acceleration (fast- forward), the 
initial state is stationary and it becomes back to the stationary state at the end of 
the fast-forward without leaving any disturbance on the WF. For the fast-forward 
of adiabatic orbital dynamics of a charged particle, we must control the driving 
filed, but we do not have to magnify the magnetic field from that of adiabatic 
dynamics, while in the standard fast-forward, the magnification of the magnetic 
field is inevitable. As typical examples, we showed fast-forward of adiabatic wave 
packet squeezing by magnetic field and adiabatic transport in EMF. Furthermore 
we showed the fast-forward of adiabatic spin dynamics in time-dependent magnetic 
field. The distinction between the present theory and Kato-Berry's transitionless 
quantum driving was elucidated. 

K.N. is grateful for the support received through a project of the Uzbek Academy of Sci- 
ences (FA-F2-084). We thanks D. Matraslov, B. Baizakov, B. Abdullaev, M. Musakhanov, 
S. Tanimura and S. Sawada for useful discussions. 



Appendix A. Driving scalar potential in systems with scalar 

potential Vq 

Here the driving scalar potential of the fast-forward of the systems with potential 
Vq is shown. The driving vector potential is given by the same form as in the case 
without Vq. 

For the fast-forward of regular (non-adiabatic) dynamics with the Hamiltonian 
Hq = 2^~(P + f Ao) 2 + Vb, we have the driving scalar potential Vff as 

VW(x,t) = a(t)V&(x,A(t)) + V$(x,t), (Al) 

where Vpp is defined by Vff in equation (|2.23[) . The driving vector potential is 
given by the same form as in equation (|2.20p . 
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The driving scalar potential for the fast-forward of adiabatic dynamics with 
Hamiltonian 



is given by 



g0= (P + Ao(x,fl(t)))» 



V FF (x,t) = Vo(x, A(t)) + V$!(x,i), (A3) 



where Vp^ is defined by Vff in equation (|3.26l) . 

In section [H we showed the fast-forward of adiabatic transport of wave packet 
trapped by magnetic field. Here we show the driving scalar potential for adiabatic 
transport of wave packet trapped by electric and magnetic field with Vo- In adiabatic 
dynamics, the trapping scalar potential is also shifted as 

V = U(x-et,y,z), (A 4) 

as well as vector potential, where U is a static trapping scalar potential. In such 
case, the driving scalar potential is represented as 

V FF (x,t) = U{x-R(A(t)),y,z) + vj$(x,t), (A 5) 

where Vpp is defined by V FF in equation (|4.17p . 

Appendix B. Fast-forward of adiabatic spin dynamics with 

adiabatic phase 

We compare the driving field in equation (I5.22p with the one obtained by Berry 
(Berry 2009) with the use of Kato's formalism (Kato 1950). Under the magnetic 
field 

sin 9 cos ip 

B (R(t)) = B | sinflsinyj ), (B 1) 

cos 9 

ci and C2 in equation (|5.11l) are represented as 

Q 

ci = cos — (B 2a) 

n 

c 2 = e^sin-, (B2b) 

respectively. Here all of Bo, 9 and tp are dependent on R(t) which is slowly chang- 
ing in time. Substituting equation (|B 2\ into equation (|5.22p . our driving field is 
expressed as 

Bf F \t) = BM(R) -2hea(~sm<p + -^<psm^ cos ^ cos <p), (B3a) 
B\f F \t) = B^(R) + 2hsa(^ R d -cos^-^sin d -cos 9 - S m^, (B3b) 
Bf F \t) = Bi"HR) + 2hea^ S ^ d -. (B3c) 
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So far, in our main text, we derived the driving field which accelerates the 
adiabatic dynamics except for the uniform phase. Now, let us choose the regularized 
standard state with the adiabatic phase as 

¥ re9 \R(t)) = ( Cl f5SJ ) e"*/o W)>«* e *W (B4) 

where £(£) is the adiabatic phase given by 



' (B5) 

*f««m+*T&- (B6) 



In this case, hn and h\2 in the regularized Hamiltonian corresponding to equation 
(|5.14D are represented as 

hu = Mcl^+c 2 ^) + L(\ Cl \ 2 -\c 2 \"), (B7a) 

hi2 = i7l( C ;||- Cl ||)+2L Cl 4 ) (B7b) 

where 

L = - lh ^m +c *m ) - (B8) 

In the analogous way used in Section [5l the driving Hamiltonian is represented as 
equation (jOTj) . With the use of equations (pT2T|) . (|B2|) . (|B6|) and (|B7)) . we obtain 
the driving field as 

Bf F \t) = B { X 0) (R) - hea(t)(— simp + sin 6> cos 6» cos <p) , (B9a) 

f)9 f)(n 

Bi FF Ht) = B^(R)+hsa( — — cos 1^—7- — sin 6 cos 9 sin if) , (B9b) 
" J oR oR 

Bi FF \t) = B^(R) + hea^sm 2 9. (B9c) 

This kind of driving field was already obtained by Berry (Berry 2009), but is 
different from ours which accelerates the adiabatic dynamics without the adiabatic 
phase. In fact, we see the polar angle 9 in equation (IB 9[) wherever 9/2 appears 
in equation (|B 3p . Interestingly, in our examples of spin inversion and LZ model 
where no adiabatic phase appears (:^| = 0), both equations (|B 3[) and (|B 9(1 lead 
to the identical driving field. In particular the y-component of Bfp in Landau- 
Majorana-Zener model (Landau 1932; Majorana 1932; Zener 1932) is available from 
both equations (|B 3|) and (|B 9ft . The distinction between equations (|B 3[) and (|B 9p 
will manifest itself in the case of the winding Landau-Zener model (Berry 1990; 
Nakamura & Rice 1994; Bouwmeester et.al. 1996), where the adiabatic phase is 
non- vanishing. 
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